33 research outputs found
Parallel solution methods and preconditioners for evolution equations
The recent development of the high performance computer platforms shows a clear trend towards heterogeneity and hierarchy. In order to utilize the computational power, particular attention must be paid to ļ¬nding new algorithms or adjust existing ones so that they better match the HPC computer architecture. In this work we consider an alternative to classical time-stepping methods based on use of time-harmonic properties and discuss solution approaches that allow eļ¬cient utilization of modern HPC resources. The method in focus is based on a truncated Fourier expansion of the solution of an evolutionary problem. The analysis is done for linear equations and it is remarked on the possibility to use two- or multilevel mesh methods for nonlinear problems, which can enable further, even higher degree of parallelization. The arising block matrix system to be solved admits a two-by-two block form with square blocks, for which a very eļ¬cient preconditioner exists. It leads to tight eigenvalue bounds for the preconditioned matrix and, hence, to a very fast convergence of a preconditioned Krylov subspace or iterative reļ¬nement method. The analytical background is shown as well as some illustrating numerical examples
Preconditioning of discrete state- and control-constrained optimal control convection-diffusion problems
We consider the iterative solution of algebraic systems, arising in optimal control problems constrained by a partial differential equation with additional box constraints on the state and the control variables, and sparsity imposed on the control. A nonsymmetric two-by-two block preconditioner is analysed and tested for a wide range of problem, regularization and discretization parameters. The constraint equation characterizes convection-diffusion processes